Fermion scattering off electroweak phase transition kink walls with hypermagnetic 

fields 
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We study the scattering of fermions off a finite width kink wall during the electroweak phase 
transition in the presence of a background hypermagnetic field. We derive and solve the Dirac 
equation for such fermions and compute the reflection and transmission coefficients for the case 
when the fermions move from the symmetric to the broken symmetry phase. We show that the 
chiral nature of the fermion coupling with the background field in the symmetric phase generates an 
axial asymmetry in the scattering processes. We discuss possible implications of such axial charge 
segregation for baryon number generation. 
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I. INTRODUCTION 

The possible existence of magnetic fields in the early 
universe has recently become the subject of intense re- 
search due to the many interesting cosmological implica- 
tions that these entail |IJ]. For instance, magnetic fields 
can influence big bang nucleosynthesis (BBN), affecting 
the primordial abundance of light elements and the rate 
of expansion of the universe. The success of the standard 
BBN scenario can be used to set limits on the strength 
of the magnetic fields at this epoch. Moreover, at decou- 
pling, long range magnetic fields can induce anisotropics 
in the cosmic microwave background radiation. Tem- 
perature anisotropies from COBE results place an upper 
bound B ~ 1 0~ 9 G for homogeneous fields (B refers to 
the intensity that the field would have today under the 
assumption of adiabatic decay due to the Hubble expan- 
sion) |2|] . In the case of inhomogeneous fields their effect 
must be searched for in the Doppler peaks Q and in the 
polarization of the CMBR @. The future CMBR satel- 
lite missions MAP and PLANCK may reach the required 
sensitivity for the detection of these last signals. 

Another interesting cosmological consequence is the ef- 
fect that primordial magnetic fields could have had on the 
dynamics of the electroweak phase transition (EWPT) at 
temperatures of the order of T ~ 100 GeV. In fact, it has 
been recently pointed out that, provided enough CP vi- 
olation exists, large scale primordial magnetic fields can 
be responsible for a stronger first order EWPT §1 [J 
(see however Ref. M ) . The situation is similar to a type 
I superconductor where the presence of an external mag- 
netic field modifies the nature of the superconducting 
phase transition due to the Meissner effect. 

Recall that for temperatures above the EWPT, the 
SU(2)xU(l)y symmetry is restored and the propagat- 
ing, non-screened vector modes that represent a mag- 
netic field correspond to the U(l)y group instead of to 
the U(l)em group, and are therefore properly called hy- 
permagnetic fields. 



In a previous work |Q| , we have shown, by using a sim- 
plified picture of a first order EWPT, that the presence 
of such fields also provides a mechanism, working in the 
same manner as the existence of additional CP violation 
within the SM, to produce an axial charge segregation in 
the scattering of fermions off the true vacuum bubbles 
The asymmetry in the scattering of fermion axial modes 
is a consequence of the chiral nature of the fermion cou- 
pling to hypermagnetic fields in the symmetric phase. 
The simplification consisted in considering the limit of 
an infinitely thin bubble wall. This assumption allowed 
us to formulate the problem in terms of solving the Dirac 
equation with a position dependent fermion mass, pro- 
portional to a step function, this last being zero in the 
false phase and non-vanishing in the broken symmetry 
phase. This treatment rendered analytic solutions from 
where reflection and transmission coefficients for axial 
modes were straightforward computed. 

In spite of the relative ease for the computation in 
such a scheme, there are two limitations related to the 
sudden change in the Higgs field profile that needed to 
be addressed. First, it is well known that the negative 
energy solutions of the Dirac equation become important 
in situations where the potential energy term changes 
over distances smaller than the particle's Compton wave 
length. Second, the height and width of the wall are 
typically related to each other in such a way that it is 
not entirely realistic to vary one without affecting the 
other. 

In this paper we overcome the above limitations by al- 
lowing a finite width of the Higgs field profile. Working 
in the thin wall regime, we use the kink solution of the 
Higgs field to formulate and solve the Dirac equation in 
the presence of hypermagnetic fields. We compute explic- 
itly transmission and reflection coefficients for the axial 
modes incident on the wall from the symmetric phase. 
Since these are related to the corresponding coefficients 
for fermions incident from the broken symmetry phase by 
CPT and Unitarity, we find that the axial charge segrega- 
tion still happens during fermion scattering of this wall. 
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The existence of such asymmetric reflection for the axial 
modes provides a bias for baryon over antibaryon produc- 
tion. In the absence of hypermagnetic fields, this mecha- 
nism has been proposed and studied in Refs. [if], [T^] 
in extensions of the SM. 

The outline of this work is as follows: In Sect. ||, we 
briefly review how the kink solution for the spatial pro- 
file of the Higgs field is obtained from a finite tempera- 
ture effective potential. In Sect. Ill, we set up the Dirac 



equation for fermions moving in this background Higgs 
field in the presence of an external hypermagnetic field. 
Section. |y| is devoted to a rather technical discussion 
about the solutions of this equation and their proper- 
ties. In Sect. 0, we use the above solutions to compute 
reflection and transmission coefficients for axial fcrmion 
modes moving from the symmetric phase toward the bro- 
ken symmetry phase. We show that these coefficients dif- 
fer for the two distinct helicity modes. Finally in Sect. |Vl, 
we conclude by looking out at the possible implications 
of such axially asymmetric fcrmion reflection and trans- 
mission. 



II. KINK SOLUTION 

To describe the EWPT, we start by writing the effec- 
tive, finite temperature Higgs potential which, including 
all the one-loop effects and ring diagrams, looks like 

A , A 



^eff ih,T) = l{T-'-T^) K l — 8 T h A 



(1) 



2 V 4 

where h = ^(WH) 1 / 2 is the strength of the SU(2) 
Higgs doublet H whose vacuum expectation value is 
given by 

v 



(H) = 



V2 



(2) 



The parameters 7, 5 and A have been computed pertur- 
batively to one loop and can be expressed in terms of v, 
the SU(2) gauge boson masses and the top mass. Their 
explicit expressions can be found elsewhere (see for exam- 
ple Ref. ||). S is the parameter responsible for the first 
order nature of the phase transition. It is the parame- 
ter that gets enhanced in the presence of hypermagnetic 
fields. T c is the critical temperature at which spinodal 
decomposition proceeds. 

We can write the effective potential in a more trans- 
parent form |l3| by introducing the dimensionless tem- 
perature d and the dimensionless Higgs field strength tp 
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At 

S T 



1-1^ 
T 



(3) 



in terms of which, the effective potential, Eq. (|l|), be- 
comes 



For simplicity, we work in the approximation where the 
energy densities of both the unbroken and broken phases 
are degenerate. This happens for a value of — 2. In 
this approximation, the phase transition is described by 
a onc-dimcnsional solution for the Higgs field, called the 
kink, which separates the two phases. This is given by 



tp(x) = 1 + tanh(x) , 
where the dimensionless position coordinate x is 

S T 



2A 



(5) 



(6) 



The parameter v2A/(d T) represents the width of the 
domain wall jlij . It can also be checked that this param- 
eter becomes smaller in the presence of hypermagnetic 
fields. 

In terms of the kink solution we can see that x — — 00 
represents the region outside the bubble, that is the re- 
gion in the symmetric phase. Conversely, for x = +00, 
the system is inside the bubble, that is in the broken 
phase. The kink wall propagates with a velocity deter- 
mined by its interactions with the surrounding plasma. 
This velocity can be anywhere between 0.1-0.9 the speed 

of light pa. 



III. DIRAC EQUATION FOR AXIAL 
FERMIONS IN A BACKGROUND 
HYPERMAGNETIC FIELD 

In the presence of an external magnetic field, we need 
to consider that fermion modes couple differently to the 
field in the broken symmetry and the symmetry restored 
phases. 

For z < 0, the coupling is chiral. Let 

*ij = ^(1 + 75)* 

*l = ^(1-75)* (7) 

represent, as usual, the right and left-handed chirality 
modes for the spinor 'J, respectively. Then, the equa- 
tions of motion for these modes, as derived from the clcc- 
troweak interaction Lagrangian, are 

(iff - ^g'4)^L-m{z)^ R = 

(i? - ^-g'4)y R -m{z)* L - 0, (8) 

where vr^l are the right and left-handed hypercharges 
corresponding to the given fermion, respectively, g' the 
U (l)y coupling constant and we take = (0, A) repre- 
senting a, not as yet specified, four-vector potential hav- 
ing non-zero components only for its spatial part, in the 
rest frame of the wall. 
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The set of Eqs. (ra) can be written as a single equation 
for the spinor ^ = ^>/j + ^ by adding up the former 
equations 



{v 



-j9 (l + 75j + — 9 (I-75) 



i{z)) 



* = 



(9) 



where the fermion mass m(z) is proportional to the vac- 
uum expectation value of the Higgs field. Hereafter, we 
explicitly work in the chiral representation of the gamma 
matrices where 



7° = 



-I 

1 



7 = 



<T 

-a 



75 = 



I 
-I 



.(10) 



Within this representation, we can write Eq. (j^j as 

{1$ - gA^ -m(z)}* = 0, (11) 
where we have introduced the matrix 



Q = 



ifg'I 

m g 'i 



(12) 



We now look at the corresponding equation in the bro- 
ken symmetry phase. For z > the coupling of the 
fermion with the external field is through the electric 
charge e and thus, the equation of motion is simply the 
Dirac equation describing an electrically charged fermion 
in a background magnetic field, namely, 



i{z)) 



* = 0. 



(13) 



In the following section, we explicitly construct the so- 
lutions to Eqs. ([ll]) and (|ll]) with a constant magnetic 
field, requiring that these match at the interface z = 0. 



IV. SOLVING THE DIRAC EQUATION 



Let us first find the solution to Eq. (|l l|) , namely, for 
fcrmions moving in the symmetric phase, z < 0. For this 
purpose, we look for a solution of the form 



(14) 



Inserting this expression into Eq. (11), we obtain 

{ - d 2 - igd'Ap - l^GF^ - 

2^ QA^ + Q 2 A^ + vfd,jn{z) } $ = , (15) 
where, as usual, 



F„ 



da A v — d v A, 



(16) 



Given that the present analysis is carried out in the thin 
wall approximation, which in physical terms means that 
the spatial region over which the fermion mass changes 
is small compared to other relevant length scales such as 
the particle mean free path, our choice of external gauge 
fields and thus field strengths should try to capture this 
information, namely, that the change in the magnetic 
field strength occurs over a small spatial region. An ex- 
act treatment of the gauge fields should be to find the 
configuration that incorporates the boundary conditions 
imposed by the change in the Higgs profile. This treat- 
ment will render continuous gauge fields across the in- 
terface. An additional feature will be the generation of a 
component of the magnetic field directed along the trans- 
verse direction, since the longitudinal component of the 
gauge fields will vary along the longitudinal direction. 
This component will be localized also in the small region 
comprised by the phase boundary. 

Notice then that when including this transverse com- 
ponent of the magnetic field into the Dirac equation, the 
separation of variables that factorizes longitudinal and 
transverse motion will not be possible. However, as long 
as this transverse field is confined to a small region and 
its strength is not too large to avoid capturing low energy 
fermions, it is a reasonable approximation to consider a 
constant field in the longitudinal direction since when the 
incident flux is not lost in modes captured on the wall, 
the probabilities of transmission and reflection depend on 
the particle currents computed in the asymptotic regions 
which in turn depend only on the fermion coupling to the 
external field, already constant in these regions. 

To estimate the magnitude of such a transverse field 
able to capture fermions on the wall, let us consider the 
following classical argument. Equilibrium between the 
Lorentz force and the centrifugal force gives for the radius 
of the orbit for a particle trapped in the wall R = p/(eB) 
where p is the particle's momentum, e is its charge and B 
the strength of the magnetic field. Taking R as the wall 
width A, and since [see Eq. @], B = b/(\ 2 ), then b = 
pX/e. Taking p ~ T (the momentum of a typical particle 
in thermal equilibrium) and since A ~ [see Eq. (](])] 
means that 6^3 which is already a very high value of 
the magnetic field. Particles with smaller momenta could 
however be trapped in this transverse component of the 
field. For the purposes of the present work, we postpone 
the consequences of such external field configuration and 
consider a piece wise constant magnetic field. 

For definiteness, let us consider the field B = Bz point- 
ing along the z direction. In this case, the vector poten- 
tial A can only have components perpendicular to z and 
the solution to Eq. (|l5|) factorizes as ll(| |l^] 



$(t,x) = T(x,y)$(t,z) 



(17) 



We concentrate on the solution describing the motion of 
positive energy fermions perpendicular to the wall, i.e., 
along the z axis. We thus look for stationary states, 
namely 



$(t,z) 



-iEt 



*(*) 



(18) 
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Therefore, working in the Lorentz gauge, d^A^ = 0, 
Eq. ( filf ) becomes 

+^ 3 ^^ + E 2 + ^BG 1 1 1 2 }^(z)=0. (19) 

Notice that Eqs. @ and have the appropriate limit 
when yn — = e, corresponding to the description of 
fcrmions coupled with their electric charge to a back- 
ground magnetic field . 

We now expand $>(z) in terms of the eigen-spinors u± 
(,s = l,2)of 7 3 @, 



(20) 



( M 













±1 


4 = 


j) 


V o) 







These spinors have the properties 



7 3 4 2 


= ±lU_t 


7°4 


= T«4 


7°4 


= ±iw^- 


7V4 


= — iu± 


12 2 

7T«± 


= +iu 2 ± 


1,2 

75M± 


1,2 
= % . 



(21) 



Writing 

$(z) = <^( Z K + + 4>\{z)u 



(j) 2 _{z)u 2 



(22) 



and inserting this expression into Eq. (|l9|), we obtain 



dz 2 



^ + (E 2 -m 2 )+g>B i - y ±± y ^ 



dm 

dz 



g'B iyL - yR) ^_(z) = 



+ -_l + { E 2 - m 2 ) +g 'B {yL+VR) 



dm 



dz 2 dz 



(23) 



and 



d dm 

dz 2 dz 



d 2 dm 
dz 2 dz 



+ (E 2 -m 2 )-g 



, {Vl + Vr) 



B 



, (Vl - VR, 



Bcb 2 _(z) = 



(E 2 - m 2 ) - g 



, (vl + Vr) 



B 



= 0. 

(24) 



Equations ( p3| ) and (|2J), represent, each, a set of two cou- 
pled second-order differential equations. The second set 
is obtained from the first one by changing B to — B. Con- 
sequently, Eqs. (^H) and the corresponding functions and 
spinors with s = 1 describe the motion of the spin compo- 
nents parallel to to magnetic field whereas Eqs. (|2~i| ) and 
the functions and spinors with s = 2, describe the mo- 
tion of the spin components antiparallel to the magnetic 
field. Notice that in the limit when yn = yh — e, each set 
of equations decouple as is the case when describing the 
interaction of fermions with the magnetic field through 
their electric charge. 

From now on, we focus on the set of Eqs. (|23|), since, 
as we have pointed out, the solutions to Eqs. ([m]) are 
obtained from those to Eqs. ( p3|) by changing B to — B. 

We now extract the dimensions writing the equations 
in terms of the dimensionless position coordinate x given 
in Eq. (^|) and, furthermore, writing them in terms of the 
new variable 



u 

obtaining, respectively, 
r d 2 1 - 2u d 



1 — tanh(a;) 



(25) 



du 2 u(l — u) du u(l — u) 4u 2 (l — u) 2 

_ £ 2 | g'b{y L + yR) ~ 
u 2 16u 2 (l — u) 2 

, g'KvL - vr) a i , v n 



l-2u d 



16u 2 (l - uf 



du 2 it(l — u) du u(l — u) 4u 2 (l — u) 2 

C 2 , g'KvL + VR)' 



u 2 16w 2 (l — u) 2 



g'KyL - vr) ,1 , v n 



(26) 



where the parameters b and £ are related to the magnetic 
field strength and the fermion mass by 



b = B 



ST 



m(x) = 



ST 



(27) 



and furthermore, £ represents twice the ratio of the 
fermion mass to the Higgs mass, e is the energy pa- 
rameter given by 

To further simplify Eqs. (prf), we try the ansatz 

4> 1 ± (u)=u a ^l-u)^ X ± (29) 
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and examine the behavior of the resulting differential 
equations near the singular points u — and u = 1. 
Assuming that the functions x± vary slowly near these 
singularities, we obtain the conditions 



Oi\ = 



e 2 + £by R - 4£ 2 = of 



k = 



U 2 + ^by L = ft 



(30) 



(31) 



Thus, to each pair of parameters, namely (a^ft) and 
(af ,ft) 1 corresponds a pair of coupled differential equa- 
tions which we write as 



' du 2 



«(l-«)^ + [c 1L -(l + 4 L +^)«]^ 



d 
du 



-a 1L b 1L 



X l L = ±f^(u) 



XL/ 



du 2 



u(l - u)-^ + [c 1R - (1 + a 1R + b 1R )u] 



d_ 

du 



i 1R b 1R 



f 1H (u) , (32) 



where we have introduced the definitions for the functions 
f 1L and f 1R given by 



/ 



XL 



c 



XR 



u(l — u) 

c 



u(l — u) 



XR , X X_ R) 



(33) 



together with that for the parameters 
cr > given by 



XL,R ,XL.R , 

Ot and 



XL.R 

l T 

<XL,R 

T 

AL,R 



L,R 



L.R 



Pi 



L,R 



■Pi 



L,R 



2a[' R + 1 . 



and the parameter £ is given by 

, = g' Avl-vr) 

^ 2 8 



(34) 



(35) 



The consistency of both sets of Eqs. ( p2[ ) requires that in 
the limit when u — ► 1, (z — ► -co), 



= -A R X 1R . 



(36) 



symmetric phase is found by looking for the scattering 
states. For our purposes, these correspond to fermions 
incident toward and reflected from the wall. There are 
two types of such solutions; those coupled with y^ and 
those coupled with y^. For an incident wave coupled 
with (yjj), the fact that the differential equations mix 
up the solutions means that the reflected wave will also 
include a component coupled with yn (j/l). Let us clas- 
sify the solutions according to the type of wave that is 
incident toward the wall. For an incident wave coupled 
with yL, which we call type (a), the most general solu- 
tions </>±°' (u) can be written as 



4 (a V) = (4>i L y + (Af:)(<j>¥)"±(A?)(<t>i»)" 



r 



b lR 



(37) 



whereas for an incident wave coupled with yn, which we 
call type (6), the most general solutions cjy^iu) is writ- 
ten as 



& W («) = ±(^)'±(Af)(^)" 



XL\II 



b iRy 



(38) 



where the functions ((p^' 11 ) 1 and (<j) 1 ± ,R ) 11 are the lin- 
early independent solutions corresponding to the solu- 
tions (x]^' 11 ) 1 and (x]^' 11 ) 11 of the homogeneous hiper- 
geometric differential equation expressed as expansions 
around u — 1, as appropriate for the symmetric phase, 



b XL,Ry 



(1 



tp i XL,R iXL.R 

2 Fi(a T ,b T , 

XL,R . , 



£L,Ryi 



XL.R . i 

a T ■ + b T 
2F 1 (c 1L > R - 



c 1L > R ;l 



u) 



h XL,R c lL,R. 



XL.R. 



c ^,R_ a lL,R_ b XL,R + 1]l _ u) ^ 



,1L,R\ 



W± ) part are the solutions corresponding to the par- 
ticular solutions (x^' ) part f the inhomogeneous equa- 
tions. The roman superindices I and II indicate the 
asymptotic behavior of the solution. I corresponds to an 
incoming wave (traveling to the right) whereas II cor- 
respond to a reflected wave (traveling to the left) in the 
symmetric phase, given the asymptotic behavior of the 
term (1 — when u — ► 1 such that 



(1 - u)^' 



(40) 



The particular solutions are expressed in terms of the 
functions f 1L and f 1R by the method of variation of pa- 
rameters, and their explicit expressions are 



with A^' R constants. 



To solve Eqs. (|32|), we first notice that they are inhomo- 
geneous hypergeometric differential equations. The solu- 
tion appropriate to describe the motion of fermions in the 



2a? 



(^ fl (*))WM) JJ - (<#*00) J WM) J 

s~ a i (1 - s)-P? 
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xf 1H (s)ds 



i 

2af 



xf 1L (s)ds. 



(41) 



To determine the solutions, we need knowledge of the 
functions f 1LR , which in turn are given self-consistently 
by Eqs. (|33|). By substituting the formal solutions 
X± a ' b \u) into Eqs. (33), this self-consistency is expressed 



in terms of integral equations satisfied by f 1LR , given 
explicitly by 



f 1L ' R (u) = p^\u) 



— „ 1L < R t 



+ 



1 R,L , 

2a 1 ' Ji 



K 1 



(u,s)f^( S )ds,(42) 



where we have introduced the functions p 1R < L given by 



c 



p 1L (u) = 



u(l — u) 

' {[(xl L ) 1 - (x 1 ^) 1 ] 



u(l — u) 

+ [<x?) n - bPY*] a*} 

and K 1R L given explicitly by 

1 



(43) 



K LH (u s) = 

K 1L (u,s) = — -. ^- = -r 

x " (*+(«))'(*+(*))" 

+ (^(u)) w (^(«)) J - (^>))'(^( S )) IJ 1 . 



(44) 



The solution to Eqs. (M2) is found numerically jjl 



We now turn to finding the solution to Eq. (13) 



namely, for fermions moving in the broken symmetry 
phase, z > 0. This time, we look for a solution of the 
form 



* = 



(45) 



Let us continue looking only at solutions type 1. By a 
procedure similar to that leading to Eqs. (26), the corre- 
sponding expressions for the functions <j>±(z) : represent- 
ing a transmitted wave moving to the right in the broken 
symmetry phase become [fill 



44 (u) = Bi±tt ai (l-«) A 

x 2 F 1 {a 1 T + l-c 1 ,b 1 T + l-c 1 ,2-c 1 ;u), 



(46) 



with B\± constants and where the hypergeometric func- 
tion 2-Fi is expressed as an expansion around u = 0, as is 
appropriate for this region. The parameters aL, b\- and 
c 1 are given by 



bl = 



ai + 0i + i 
ai + Pi + - 



= 2ot a + 1 



(47) 



with 



ai = -^e 2 +g>b-4£ 2 



ft = ^V^+Vb 



(48) 



Notice that since in the broken symmetry phase there 
should not be a propagating component corresponding to 
the Z° field || [|, the parameter representing the mag- 
netic field strength b' is related to b and Weinberg's angle 
9 W by 



b' 



cos 6 W 



(49) 



which in turn implies that the coupling of the fermion 
with the magnetic field is given by 



eb' = g'b . 



(50) 



The complete solution to the problem is found by match- 
ing the functions <j>± a ^ (u) and </>±(u) as well as their 
derivatives across the interface at u = 1/2. These con- 
ditions represent four algebraic complex equations that 



determine the four complex constants A^' R and B\± 



V. TRANSMISSION AND REFLECTION 
COEFFICIENTS 

The fact that the amplitudes for the axial modes in 
the symmetric phase, Eqs. (|37| ) and (|38| ) are not the 
same, means that there is the possibility of building an 
axial asymmetry during the scattering of fermions off 
the wall. To quantify the asymmetry, we need to com- 
pute the corresponding reflection and transmission coef- 
ficients. These are built from the reflected, transmitted 
and incident currents of each type. Recall that for a given 
spinor wave function \& ', the current normal to the wall 
is given by 



j = ^tyy^ 



(51) 



The currents need be computed in the asymptotic regions 
far away from the wall where the amplitudes represent 
plane waves with well defined direction of motion. 

We now prepare the incident fermion from the symmet- 
ric phase in such a way that when coupled with a given 
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chirality (left-handed for waves type (a), right-handed for 
waves type (6)), it corresponds to the same helicity. Since 
in the symmetric phase the fermion mass is asymptoti- 
cally zero, both the chirality and helicity operators can 
be simultaneously defined and their eigenvalues coincide. 
This is no longer the case when the fermion moves in the 
broken symmetry phase where its mass is different from 
zero. Nevertheless, since scattering off the wall does not 
change the direction of the fermion spin (modes 1 and 
2 evolve independently) the fermion helicity is preserved 
during transmission and reversed upon reflection. 

For left-handed incoming waves, the incident current 
Jj nc (lower case indexes I and r denote helicity modes) 
is thus given by 

4 lc = 4|ze + 2/3 2 L | 2 (52) 

whereas the reflected and transmitted currents J r ' ~, J| ra 
are given respectively by 

J r Tei = 'l{\A5\ 2 \-ie + 2rt\ 2 -\AK\ 2 \te + 2[3?\ 2 } 

J{ ra = {\B 2+ [2(£ - a 2 ) - ie] - £ 2 _[2(£ + a 2 ) + ie]\ 2 

+ |£ 2 +[2(£ - a 2 ) + ie] - B 2 -[2(£ + a 2 ) - ie]| 2 } . 

(53) 

On the other hand, for right-handed incoming waves, the 
incident current JjL„is thus given by 

Jf nc = 4|ze + 2/3f| 2 (54) 

and the reflected and transmitted currents j, ^ ra are 
given respectively by 

4f - 4 {Kl 2 H e + 2 #f - K I 2 1« + 2/3f| 2 } 

J t r ra = {\B 1+ [2(£ - ai) - ie] + B x _[2(£ + a x ) + ie]\ 2 

- \B 1+ [2{i - ai ) + ie] + Bi_[2(£ + a t ) - ze]| 2 } . 

(55) 

For a left-handed incident particle, the reflection and 
transmission coefficients are given as the ratios of the 
corresponding reflected and transmitted currents, to the 
incident one, respectively, projected along a unit vector 
normal to the wall. These are 

Rl^r = -^ref/ J inc 

Tl-M = 4a/4c ( 56 ) 

The corresponding coefficients for the axially conjugate 
process are 

Rr^i = - J leil J lnc 

T r ^ r = JtraMnc ■ ( 57 ) 

Figure 1 shows the coefficients R^ r and R r ^i as a func- 
tion of the magnetic field parameter b for a value of twice 
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FIG. 1: Coefficients Ri^ r and R r ^i as a function of the mag- 
netic field parameter b for £ = 3.5, e — 7.03, yn = 4/3, 
— 1/3. The value for the U(l)y coupling constant is taken 
as g' — 0.344, corresponding to the EWPT epoch. The dots 
represent the computed values. 

the ratio of fermion to Higgs mass £ = 3.5, an energy pa- 
rameter e = 7.03, hypercharge values ya — 4/3, yx, = 1/3 
and for a value of g' — 0.344, as appropriate for the 
EWPT epoch. Notice that when 6^0, these coeffi- 
cients approach each other and that the difference grows 
with increasing field strength. 

Figure 2 shows the reflection and transmission coeffi- 
cients as a function of the energy parameter e scaled by 
twice the height of the barrier 2£. Figure 2a shows the co- 
efficients Ri^r and Ti^i and Fig. 2b the coefficients i? r _>; 
and T r _> r for b = 0.5 and £ = 3.5, yn — 4/3, Ul = 1/3, 
g' = 0.344. Since the solutions in Eqs. ( (46[) are com- 
puted assuming that the transmitted waves are not ex- 
ponentially damped, their energies have to be taken such 
that the parameters a x ,2 are imaginary which in turn im- 
plies that for waves type 1, e > \J 4£ 2 — g'b whereas for 
waves type 2, e > -y/4£ 2 + g'b. It can be checked that 
R r ^i + T r ^ r = 1 and Ri^ r +T;_ > ; = 1 within the numer- 
ical precision of the calculation, which means that the 
analysis respects unitarity. 



VI. CONCLUSIONS 

In this paper we have derived and solved the Dirac 
equation for fermions scattering off a first order EWPT 
bubble wall with a finite width in the presence of a mag- 
netic field directed along the fermion direction of motion. 
In the symmetric phase, the fermions couple chirally to 
the magnetic field, which receives the name of hypermag- 
netic, given that it belongs to the U(1)y group. We have 
shown that the chiral nature of this coupling implies that 
it is possible to build an axial asymmetry during the scat- 
tering of fermions off the wall. We have computed re- 
flection and transmission coefficients showing explicitly 
that they differ for left and right-handed incident parti- 
cles from the symmetric phase. The results of this more 
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FIG. 2: Reflection and transmission coefficients as a function 
of the energy parameter e scaled by twice the height of the 
barrier 2£ for b — 0.5 and £ = 3.5, yn = 4/3, i/l ~ 1/3, 
g' = 0.344. Figure 2a (upper panel) shows the coefficients for 
incident, left-handed helicity modes and Fig. 2b (lower panel) 
for incident, right-handed helicity modes. In both figures, the 
dots represent the computed values. 



realistic, albeit numerical calculation where we allow for 
a finite wall width are in qualitative and quantitative 
agreement with those previously found in Ref . |9| , where 
the wall was modeled as a step function. 

It could be thought that the asymmetric reflection 
found in this work could be washed out when considering 
the averaging over the different angles of incidence of the 
fcrmion flux. This is not the case as we proceed to show. 
Let us first look a the situation in which the direction 
of the magnetic field is reversed with respect to the case 
studied here. This a physically relevant scenario since 
during the phase transition, fermions are scattered on 
opposite sides of the bubbles and if the sign of the asym- 
metry would depend on the direction of the magnetic 
field with respect to the direction of fermion incidence, 
then the building of an axial charge density in one side of 
the bubble would compensate the building of this charge 
on the other side, thereby canceling the effect. However, 
it is easy to convince oneself that this is not the case. By 



looking at Eqs. ( p3| ) and (|24|), we see that changing B 
to — B interchanges one set of equations with the other, 
leaving intact the coupling. Physically this is also easy to 
understand since the fermion coupling with the external 
field is through its spin. Changing the direction of the 
field exchanges the role of each spin component but since 
each chirality mode contains both spin orientations, it 
does not affect the final probabilities. 

Now suppose that the original direction of motion of 
the fermion is not parallel to the direction of the mag- 
netic field and therefore its velocity vector contains a 
component perpendicular to the direction of the field. In 
this case, due to the Lorentz force, the particle circles 
around the field lines maintaining its velocity along the 
direction of the field. The motion of the particle is thus 
described as an overall displacement along the field lines 
superimposed to a circular motion around these lines. 
In the three dimensional quantum mechanical treatment 
of the problem, these circles correspond to the different 
Landau levels. We see that the originally different an- 
gles of incidence all result in the same overall direction 
of incidence. Nonetheless, it is certainly true that these 
circular trajectories could be regarded as the paths where 
the wave function of the particle picks up a phase in the 
same manner as in the Aharanov-Bohm effect. However, 
since there is no definite phase relation of the incident 
fermions, these phases have to be regarded as randomly 
distributed. Thus, the addition of the wave functions 
at the interference point (minus infinity for the reflected 
waves and plus infinity for the transmitted waves) has to 
be done incoherently which precludes any possible de- 
structive effect of these phases on the overall particle 
fluxes. 

We also emphasize that, under the very general as- 
sumptions of CPT invariance and unitarity, the total ax- 
ial asymmetry (which includes contributions both from 
particles and antiparticles) is quantified in terms of the 
particle (axial) asymmetry. Let pi represent the number 
density for species i. The net densities in left-handed 
and right-handed axial charges are obtained by taking 
the differences Pl — Pi an d PR — PR, respectively. It is 
straightforward to show that CPT invariance and 
unitarity imply that the above net densities are given by 



PL 

PR 



PL 
PR 



(f 



f b )(Rr- 



Ri- 



(58) 



where f s and f b are the statistical distributions for par- 
ticles or antiparticles (since the chemical potentials are 
assumed to be zero or small compared to the temper- 
ature, these distributions are the same for particles or 
antiparticles) in the symmetric and the broken symme- 
try phases, respectively. From Eq. (|58|), the asymmetry 
in the axial charge density is finally given by 

(PL - PL) - (PR - PR) = 2(/ S - f)(Rr^l - R^r).(59) 

This asymmetry, built on either side of the wall, is disso- 
ciated from non-conserving baryon number processes and 
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can subsequently be converted to baryon number in the 
broken symmetry phase where sphaleron induced transi- 
tions are taking place with a large rate. This mechanism 
receives the name of non-local baryogenesis |fio| , O , |l^ , |2C[| 
and, in the absence of the external field, it can only be re- 
alized in extensions of the SM where a source of CP viola- 
tion is introduced ad hoc into a complex, space-dependent 
phase of the Higgs field during the development of the 
EWPT @. 

Since another consequence of the existence of an ex- 
ternal magnetic field is the lowering of the barrier be- 
tween topologically inequivalent vacua p2[ , due to the 
sphaleron dipolc moment, the use of the mechanism dis- 
cussed in this work to possibly generate a baryon asym- 
metry is not as straightforward. Nonetheless, if such pri- 
mordial fields indeed existed during the EWPT epoch 



and the phase transition was first order, as is the case, 
for instance, in minimal extensions of the SM, the mech- 
anism advocated in this work has to be considered as 
acting in the same manner as a source of CP violation 
that can have important consequences for the generation 
of a baryon number. These matters will be the subject 
of an upcoming work ]23| . 
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